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We analytically investigate the ground-state properties of two-component Bose-Einstein conden- 
sates with few 87 Rb atoms inside a high-quality cavity quantum electrodynamics. In the SU(2) 
representation for atom, this quantum system can be realized a generalized Dicke model with a 
quadratic term arising from the interatomic interactions, which can be controlled experimentally 
by Feshbach resonance technique. Moreover, this weak interspecies interaction can give rise to an 
important zero-temperature quantum phase transition from the normal to the superradiant phases, 
where the atomic ensemble in the normal phase is collectively unexcited while is macroscopically 
excited with coherent radiations in the superradiant phase. Finally, we propose to observe this 
predicted quantum phase transition by measuring the direct and striking signatures of the photon 
field in terms of a heterodyne detector out of the cavity. 
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Cavity quantum electrodynamics has an essential mi- 
crowave photon field to explore the important micro- 
scopic quantum phenomena in quantum optics since the 
strong atom-field coupling strength, which dominates 
over the dissipative losses of the quantum system, has 
been achieved]!! H Q ■ Recently, much progress has been 
made in the observation of the motional dynamics 0, [f| 
as well as the trapping^, 0] and coolingd, Q of single 
atom within the cavity mode. This also provides an in- 
teresting way to process quantum information and im- 
plement quantum computing such as nonclassical light 
sources [13, [ll| and quantum state transfer [l^. 

Recently, Bose-Einstein condensate (BEC) coupled 
with a high-quality cavity quantum electrodynamics has 
been attracted much attentions in both theory and ex- 
periment. Despite BEC is a central goal for atom 
chips, it exhibits some novel macroscopic quantum ef- 
fects, which has never predicted and observed in tradi- 
tional atom. For example, atom-atom or atom-photon 
entanglement and squeezing can be produced in the lin- 
ear regime of collective atomic recoil lasing, where the 
ground state of the condensate remains approximately 
undepleted[l3l. [lil fl5|. By solving the three- mode mas- 
ter equation in the Wigner representation, three-mode 
entanglement as well as two-mode atom-atom and atom- 
radiation entanglement can be generally robust against 
losses and decoherence, and thus regarded as a good 
candidate for the experimental observation of entan- 
glement in condensate systems [16]. Moreover, a novel 
cavity-mediated long-range atom-atom interactions has 
been also realized theoretically [l7j- However, the ex- 
perimental process has only been considered recently 
[3 H S3, El [H, El- The main difficulties arise from 
adverse vacuum requisites and sophisticated topological 
requirements on both of these art technologies. Here 
we theoretically discuss two-component BECs interact- 
ing with a high-quality cavity quantum electrodynam- 



ics. It has been shown that in two-component BECs, 
the weak nonlinear interactions controlled experimen- 
tally by Feshbach resonance technique [24| can give rise to 
a number of novel phenomena including instab ility 25 1 , 
metastabilitv [26| . wave mixing [27l . solitonJH l29l 3C|, 
dynamical bifurcation [31j|. chaos [32l l33fl . and Heisenberg- 
limited Mach-Zahnder interferometry[34j. Very recently, 
a novel second-order quantum phase transition from the 
normal to the tunneling phases has been predicted in the 
resonant case when two-component BECs are coupled 
with a periodically driven laser field [35|. 

In this paper we mainly investigate the ground-state 
properties of two-component BECs with few 87 Rb atoms 
inside an ultrahigh finesse optical cavity, which can sup- 
port a single-mode photon. We show that in the SU(2) 
representation for atom, this quantum system can be re- 
alized a generalized Dicke model with a quadratic term 
arising from the mean-field interatomic interactions. By 
means of Holstein-Primakoff transformation [36[ and bo- 
son expansion method[37[, the ground-state energy can 
be approximately evaluated and show that the weak in- 
terspecies interaction can lead to an important zero- 
temperature quantum phase transition from the normal 
to the superradiant phases, where the atomic ensemble in 
the normal phase is collectively unexcited while is macro- 
scopically excited with coherent radiations in the super- 
radiant phase. Finally, we propose to observe this pre- 
dicted quantum phase transition by detecting the direct 
and striking signatures of the photon field such as the 
well-measured intracavity intensity / oc \(a) | 2 in terms of 
a heterodyne detector out of the cavity [38|. Moreover, 
its essential features of quantum criticality such as the 
scaling behavior, critical exponent and universality are 
also given. 

A physical realization shown in Fig.l is that two con- 
densates in different hyperfine levels \F — l,mt = — 1) 
(|1)) and \F — 2,mj — 1) (|2)) are confined in a time- 
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average, orbiting potential magnetic trap. Usually, in 
the 87 Rb system, an external laser is applied to induce a 
Josephson-like coupling and the detuning of the laser is 
adiabatically changed to produce various transitions be- 
tween |1) and 1 2). However, here this tunable laser can be 
displaced by the quantized field of the cavity mode after 
two condensates are sent into the optical cavity. For a 
high-quality cavity with the length 178 pin and the mode 
waist radius 26 pm, the finesse is given by 3 x 10 5 and the 
maximum coupling strength between the 87 Rb atom and 
the cavity field is g = 2ir x 10.4 MHz, which is larger than 
the cavity field decay rate k = 2ttx 1.4 MHz and the atom 
dipole decay rate j = 2n X 1.4 MHzH^HJ, H2] ■ Thus, in 
such a cavity the single-mode photon can be considered 
and the quantum dissipation effect can be also negligible. 
Finally, in the formalism of second quantization based 
on the self-consistent Gross-Pitaevskii equations for the 
two-component BECs, the total Hamiltonian for Fig.l 
without rotating-wave approximation can be written as 

H = H p h + H at + H at -at + H a t- P h (1) 

with 

Hph = ua + a, (2) 

H at =J2 [d 3 r{*t(r)i-^ + V l (r)}* t (r)+ (3) 

1 = 1.2 m 

|*+(r)*+(r)*,(r)*l(r)}, 

H a t^at = J d 3 r gi , 2 *+(r)tf+(r)*i(r)*2(r), (4) 

E at - Vh = »«1 J d 3 r[*+(r) exp(ik • r)* 2 (r)+ (5) 

*2 (r) exp(-ik • r)* 1 (r)](a - a + ), 

where a and a + are the annihilation and creation op- 
erators of the cavity mode with the frequency u>; ^(r) 
is the boson field operator; Vi(r) is the single magnetic 
trapped potential with the frequencies u>i(i — x,y,z); 
qi = 4irh 2 pi/m is the intraspecies interactions with p\ 
being the intraspecies s— wave scattering length and m 
being the atomic mass; q\ 2 = AttH 2 pi^/m is the in- 
terspecies interactions with p\_ 2 being the interspecies 
s— wave scattering length; Q is the atom-cavity coupling 
constant; and k is the wave vector of the quantized field 
for the cavity mode. 

It has been shown that the well-known two-mode ap- 
proximation, which is defined as 'l'i(r) = ci^>i(r) and 
^ 2 (r) = c 2 (f> 2 (r) with c\ and C2 being the annihilation 



boson operators, can be used to simplify the Hamilto- 
nian (1) as 

H = ^2 faciei + -^cfcfciq) + xctc 1 c^c 2 + (6) 
1=1,2 

X(cfc 2 + c% Ci)(a + a + ) + iva + a, 

where = f d s r{^(v)[-^A + ^(r)]0 ; (r), 

m = qi Jd 3 r\Urt, X = gi, 2 /d 3 r|^(r)| 2 |^(r)| 2 
and A = ihft J c? 3 r^(r) exp(ik • r)(f> 2 (r) = 
ihil J rf 3 r02(r) exp(— ik • r)0i(r) are the atom-cavity 
interacting constant and assumed to be real for the 
sake of simplicity. However, it should be noticed that 
this two-mode approximation can o nly be applied for 
a small number of condensed atoms [33. |40| . A simple 
estimate shows that the number of atoms N should 
satisfy the condition such that Np < Tq, where p is a 
typical scattering length and rg is a measure of the trap 
size. If we take into account large trap size r$ = 10 
pm and the typical scattering length p = 5 run [41J, the 
two- mode approximation is valid for N < 2000. 

The Hamiltonian (2) can be further simplified by ap- 
plying the Schwinger relations defined as S x = (c^C2 + 
ctci)/2, S y = {c\c 2 - Cq ci)/2«, and S z = (c+ci - 
c 2 c 2)/2, where the Casimir invariant is S 2 = N(N/2 + 
l)/2. The eigenvalues m of operator S z represent the 
difference 2(Ni — N 2 ) in the number of atoms in different 
hyperfine levels, while S x and S y take on the meaning of 
relative phase between two condensates. Besides a triv- 
ial constant, the Hamiltonian (6) can be reduced to a 
generalized Dicke Hamiltonian with a quadratic term 

H = qS 2 z + u Q S z + 2XS X (a + + a) + ua+a, (7) 

where q = [(771 + rj 2 )/2 — x] and uj = ^1 - u 2 + [N — 
1)(»?2 — r?i)/2 ~ N(i] 2 —T]i)/2 for the single trapped poten- 
tial. The Hamiltonian (7) describes the collective quan- 
tum dynamics of two-component BECs with few 87 Rb 
atoms inside an ultrahigh finesse optical cavity. It can 
be seen clearly that the nonlinear quadratic term arises 
from the mean-field interatomic interactions, which can 
be controlled by the s— wave scattering lengths of atoms 
via Feshbach resonance technique. For simplicity, here 
we assume that the weak nonlinear parameter q depends 
only on the interspecies interaction, namely, on the pa- 
rameter x (or pi, 2 ) and consider the case of loq > 0, that 
is, P2 > Pi- If the s— wave scattering lengths of atoms 
are chosen as p±. 2 = (pi + p 2 )/2, the Hamiltonian (7) 
can be mapped into the standard Dicke-like Hamiltonian 
H D = Lua+a + uj S z + 2\S x (a+ + a)^. 

In order to effectively describe the collective dynami- 
cal behavior of the Hamiltonian (7), we should evaluate 
its ground-state energy. Since here the trapped atom 
number is considered to be of order 10 3 , we can use the 
mean- field approximation to arrive at the target. Fol- 
lowing the procedure of Ref . [43| . we firstly employ the 
well-known Holstein-Primakoff transformation, which is 
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denned as S+ = b+^N - b+b, S_ = y/N - b+bb, and 
S z = (b+b - N/2) with [6,6+] = 136], to rewrite the 
Hamiltonian (7), apart from a trivial constant, as 



H = Lua+a + u {b+b- N/2) + -^={b ~ ■ ■ ( M ) 



VN -b+bb)(a+ +a) + — {b+VN - b+b){\/N - b+bb) 



where v = —Nq, g = \?N\ and ujq = ujq + q ~ wo 
since the parameter q is a negligibly small number com- 
pared with the parameter loq. Secondly, we introduce 
the new boson operators by setting c+ = a+ + y/N a and 
d+ = b+ — to describe the collective behavior of 

the Hamiltonian (8). At last, in terms of boson expan- 
sion method based on the introduced operators c+ and 
d+ , the Hamiltonian (8) can be expanded by 37 



H = NH + N 1/2 H X + N°H 2 + ■■■ 



(9) 



with Hq = ujo 2 + lu {/3 2 - 1/2) - 4gaP\/k + vP 2 yfk, 
Hi = {-Lua + 2g/3Vk)(c+ + c) + [lo /3 - (2ga - 
vpVk){Vk - (3 2 /Vk)](d+ + d), H 2 = ujc+c + (cj + 
4gal3/Vk+iyk-3vp 2 )d+d+{gaf3/Vk-ijf3 2 )[(d+) 2 +d 2 } + 
(ga(3 3 /2k^ 2 )(d++d) 2 + (gVk- /3 2 /Vk){c+ +c)(d+ +d), 
where k = 1 — (3 2 . 

The critical transition point can be derived from the 
condition Hi = 0, which means that for any quantum 
system the free energy should be fixed at the minimum 
value. Therefore, two different cases for 5 > 1 and S < 1 
with S = ujuj /N(4:X 2 + ojq) must be considered. The 
corresponding auxiliary parameters a and (3 are given by 



0, S > 1 



S < 1 



,/3 = 



0, S > 1 



(10) 



and the interesting critical transition point is given by 
(5=1, namely, 



pi)/y/2Trd x dyd z m ~ 2866 Hz. If the frequency of the 
photon and the atom-field coupling strength are consid- 
ered by w = 4 x 10 8 MHz and A = 2tt x 5 MHz[Ii|,[2l|,[2l], 
respectively, the critical transition point is evaluated as 
q c = —7 Hz, whose corresponding interspecies s— wave 
scattering length is given by (pi^)c — 7.14 nm. It seems 
that this critical parameter {p\^)c can be achieved by 
current Feshbach resonance technique, which awaits ex- 
perimental validation. 

In the case of pi,2 > [pi,2)c (5 > 1), the ground-state 
energy is given by (Eo) nor ~ —Nujo/2, which depends 
only on the parameter of the BECs. However, in the 
case of pi 2 < {pi 2 )c (<5 < 1), the ground-state energy 
becomes (E a ) sup = -N[N(\ 2 /u + g/4)(l - S 2 ) +uq5/2], 
which is dependent of the parameters of both the BECs 
and the cavity field. The physics may be understood as 
follows. For the case pi t 2 > {pi, 2)0 the quadric term 
(qS 2 ) in the Hamiltonian (7) plays an important role, 
which means that the transition between two conden- 
sates is suppressed and therefore no collective excitation 
happens. However, for the other case pi^ < (^1,2)0 
the energy of the term S x {a+ + a) in the Hamilto- 
nian (7) dominates, which implies that both the internal 
Josephson tunneling and the macroscopic collective ex- 
citation occur. Therefore, we can call the normal phase 
when pi : 2 > {pi,2)c and the superradiant phase when 
Pi, 2 < (pi,a)c< Fig. 2 shows the scaled ground-state en- 
ergy Eq/N and its second-order derivative with respect 
to the intraspecies s— wave scattering length pi.2 as a 
function of pi t %. It can be seen that the second-order 
derivative of Eq /N possesses a discontinuity at the tran- 
sition point {pi t 2) C j which clearly illustrates the nature of 
second-order phase transition (the first-order derivative 
of Eq/N is continuous at the transition point). Fig. 3 
shows the scaled ground-state atom population between 
two condensates 



~ N{£$+m) ' Pl 2 < (P^*)" (12) 
-1, pi, 2 > (pi, 2 )c 



qc= N-~- 

Since the many-body interactions in the two-mode ap- 
proximation produce only small modification of the 
ground-state properties, the wavefunctions of the macro- 
scopic condensate states for the single magnetic trap can 
be given, if the acting role of the gravity is omitted, 
by 0,(r) = 7T-^{d x d y d z )- 1 / 2 exp[-(x 2 /d 2 + y 2 / d 2 y + 
z 2 /d 2 )/2] with d x — ^/h/moj x , d y = ^/h/moj y and 
d z = y/h/muj z . When the relevant parameters are cho- 
sen as lu x = ujy = 2tt x 290 Hz, uj z = 2tt x 450 Hz, 
to = 1.45 x 10~ 25 kg, pi — 3.70 nm, p 2 — 5.70 nm, and 
N = 1000, respectively, the parameter loq in the unit of 
h can be approximately calculated as u>q ~ NH 2 {p2 — 



and its first-order derivative with respect to the in- 
traspecies s— wave scattering length pi 2 as a function 
of pi ; 2, which also demonstrates the superradiant phase 
transition. 

In the rest part of this paper we discuss how to observe 
this predicted phase transition in current experimental 
setups. In general, the many-body quantum pseudospin 
state in this quantum system is not accessible to observe 
quantum phase transition. However, here we propose 
to detect the direct and striking signatures of the pho- 
ton field such as the well-measured intracavity intensity 
/ cx I (a) 1 2 in terms of a heterodyne detector out of the 
cavity [38| . In terms of the auxiliary parameters a given 
by Eq.(10), the scaled ground-state intracavity intensity 
is evaluated as 
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JV W 2 (4A 2 +W9) 2 , Pl,2 < [Pl,2)c 

), Pi, 2 > (Pl,2)c 



(13) 



Fig. 4 shows the scaled ground-state intracavity inten- 
sity I/N and its first-order derivative with respect to 
the intraspecies s— wave scattering length pi.2 as a func- 
tion of pi.2- It is interesting that this predicted quan- 
tum phase transition characterized by the non-analyticity 
of the scaled ground-state intracavity intensity I/N is 
remarkably of the first-order. When the intraspecies 
s— wave scattering length approaches the critical 
value (pi,2)c, the scaled ground-state intracavity inten- 
sity I/N vanishes as 



N 



\pi 



(Pl,2)c] ~ \P1,2 ~ (Pl,2)c| 



(14) 



Since the diverging characteristic length scale is £ ~ 
\Pi,2 ~ (Pi,2)c| with v = 1/2, the critical exponent 
for the scaled ground-state intracavity intensity can 
be derived from I/N ~ \pi^ ~ (pi,2)c\ ZV by z = 2, 



which shows the universality principle of quantum phase 
transition [44|. 

In conclusion, we have analytically discussed the 
ground-state properties of two-component BECs of 87 Rb 
atoms inside an ultrahigh finesse optical cavity support- 
ing a single-mode photon. It has been shown that the 
weak interspecies mean-field interaction can give rise to 
a novel zero-temperature quantum phase transition from 
the normal to the superradiant phases, which means that 
the weak microscopic nonlinear interaction can lead to a 
strong behavior in macroscopic scale. Moreover, we have 
proposed to observe this predicted quantum phase tran- 
sition by measuring the direct and striking signatures of 
the photon field such as the well-measured intracavity 
intensity / oc \(a)\ 2 in terms of a heterodyne detector 
out of the cavity. It is also interesting that this scaled 
ground-state intracavity intensity I/N is remarkably of 
the first-order at the transition point (pi,2)c- 
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Figure Captions 

Fig.l (Color online) Schematic diagram for two- 
component BECs of few 87 Rb atoms inside a high-quality 
cavity quantum electrodynamics. Initially, two conden- 
sates in different hyperfme levels \F = l,m/ = —1) and 
\F = 1, mf = 1) are confined in a time-average, orbiting 
potential magnetic trap. After BECs are sent into the 



optical cavity, the quantized field of the cavity mode can 
be used to adiabatically control the various transition 
between these two condensates. The heterodyne detector 
(HD) out of the cavity can be applied to conveniently 
detect the well-measured intracavity intensity. 
Fig. 2 (Color online) The scaled ground-state energy 
Eo/N versus the intraspecies s— wave scattering length 
pi,2- The relative parameters are chosen as u> x — u> y — 
2tt' x 290 Hz, lo z = 2tt x 450 Hz, m = 1.45 x 10~ 25 kg, 
pi = 3.70 nm, p 2 = 5.70 nm, w = 4 x 10 s MHz, A = 2tt x 5 
MHz and N = 1000. The critical interspecies s— wave 
scattering length can be evaluated by (pi, 2)0 = 7.14 nm. 
Insert: the second-order derivative of Eo/N versus with 
respect to pi,2 versus p\ t i. 

Fig. 3 (Color online) The scaled ground-state atom pop- 
ulation between two condensates AN/N versus the in- 
traspecies s— wave scattering length pi, 2 with the same 
parameters as those in Fig. 2. Insert: the first-order 
derivative of AN/N with respect to pi,2 versus pi, 2- 
Fig.4 (Color online) The scaled ground-state intracavity 
intensity I/N versus the intraspecies s— wave scattering 
length pi,2 with the same parameter as those in Fig. 2. 
Insert: the first-order derivative of I/N with respect to 
Pi,2 versus pi,2- 



